Edge-induced spin polarization in two-dimensional electron gas 
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We characterize the role of the spin-orbit coupling between electrons and the confining potential of the edge 
in nonequilibrium 2D homogeneous electronic gas. We derive a simple analytical result for the magnitude of the 
current induced spin polarization at the edge and prove that it is independent of the details of the confinement 
edge potential and the electronic density within realistic values of the parameters of the considered models. 
While the amplitude of the spin accumulation is comparable to the experimental values of extrinsic spin Hall 
effect in similar samples, the spatial extent of edge induced effect is restricted to the distances of the order of 
Fermi wavelength (~ lOnm). 
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I. INTRODUCTION 

One of the exciting new discoveries in solid state physics 
in the last few years has been the experimental observation 
of the extrinsic spin Hall effect 1 " 4 in GaAs heterostmctures. 
The mechanism of this effect relies on the spin-orbit cou- 
pling between the spin of electrons with the perturbing po- 
tential of impurities. Another example of similar coupling 
is the Rashba-Bytchkov interaction in asymmetrically doped 
GaAs quantum wells containing two-dimensional electron gas 
(2DEG) where the coupled potential comes from the internal 
electrostatic electric field induced by the structural asymmetry 
of the well. In modeling both of these situations, the atomic 
potential of the ideal bulk semiconductor enters only through 
the renormalization of the effective mass and the spin-orbit 
interaction strength. 

Interestingly, the Rashba-Bytchkov interaction has been 
also suggested as a source of spin Hall-like phenomenology 5 
for a clean 2DEG. While the effect has been shown to disap- 
pear in extended 2D systems with arbitrarily weak disorder 6,7 , 
it is now known that finite size of the sample in combina- 
tion with Rashba-Bytchkov coupling results in the transverse 
spin current and accumulation of opposite spin densities at the 
edges of the sample 8,9 . Furthermore, it has been found that 
boundaries might directly influence the spin polarization due 
to the Rashba-Bytchov coupling either in very wide 2D elec- 
tronic systems with an edgol£t~— or narrow Rashba strip a 15 i 16 . 
However, the physical origin of the latter is quite different 
from the former: whereas the bulk Rashba-Bytchkov effect 
arises from accelerations of electrons in the external electric 
field, the latter result from the reflection of electrons from the 
sample's boundary. 

Motivated by this development, it is natural to explore fur- 
ther alterations of the perfectly periodic bulk potential and 
its consequences on the spin polarization in the presence of 
the current. Several authors have considered such a situation 
in quantum wires with a parabolic confining potential^ - —, 
wider strips with parabolic 21 or abrupt 22 confinement edge. 
These alterations in the potential landscape can be also viewed 
as the "impurity" which leads to nonzero spin polarization in- 
duced by the flow of the electric current. Recently, exper- 



imental evidence of this effect of the in-plane field on spin 
polarization has been reported 23 . Further enhancement of this 
kind of spin-polarization has been suggested using a transport 
through chaotic quantum dot 24 or resonant tunable scattering 



center—. 

In our previous work 2 ^ we have compared the magnitude 
of the spin polarization induced by the latter mechanism with 
the one induced by the Rashba-Bytchkov coupling and scat- 
tering of the edge 10 . It turns out that for a typical 2DEG with 
an edge, it is about three orders of magnitude larger than the 
one caused by the edge and the Rashba-Bytchkov interaction. 
This substantial difference is mostly due to the fact that the pa- 
rameter characterizing spin-orbit coupling appears in second 
order in the case of the Rashba-Bytchov-based mechanism 10 
whereas for the edge-induced spin-orbit coupling it comes al- 
ready in the first order. 

In this paper we consider the edge of the 2DEG in Si doped 
GaAs quantum well, a system for which the spin Hall effect 
has been experimentally observed. This allows us to estimate 
the values of all the parameters entering and we can compare 
the importance of this edge-effect to the contribution of other 
mechanisms leading to current-induced spin polarization. The 
models considered here are reasonably realistic yet simple 
enough to be solved almost analytically. We derive a gen- 
erally valid simple analytical formula for the induced number 
of electrons with un-compensated spin per unit length of the 
edge, 



hq, 



(D 



where j is the electrical current density in the 2DEG, Ofc is 
the parameter giving the strength of the spin-orbit coupling 26 , 
and m e and q e are the effective mass and the charge of elec- 
tron respectively. Noticeably, this result is independent of the 
density of the 2DEG and the details of the edge potential as 
long as the electrons are confined within the sample. 
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II. THE MODEL OF THE EDGE AND SPIN-ORBIT 
COUPLING 

Let us consider a sample with 2DEG created within a 
GaAs-based quantum well. The electrical current in the plane 
of the 2DEG will be driven along the y direction. We will 
be interested in the physics close to one of the two edges of 
the sample along which the current flows. The coordinate in 
the direction perpendicular to the edge will be x, x > corre- 
sponding to the region where the 2DEG is present (see Fig. [TJ. 

We will describe the electronic states for electrons in 
the 2DEG within the effective mass approximation, treat- 
ing the electrons as noninteracting quasiparticles. Qualita- 
tive changes in our results introduced within a self-consistent 
mean-field treatment will be studied afterwards. The length 
scale characterizing the electrons is given by their Fermi 
wavelength. It typically attains values 33 Xf = ^Ik/hid ~ 
2.5a* B , where «2D is the two-dimensional electronic density 
and a* B = 9.79nm is the effective Bohr radius in GaAs. Since 
this is about two orders of magnitude larger than the inter- 
atomic distances, we can employ the effective mass theory 
and approximate the form of the confining potential with a 
simple functional form. 

In our work will assume that the confining edge potential, 
V(x), is independent of the coordinate y, directed along the 
edge. Within our analytical derivations we model V(x) as an 
abrupt step of height AV, V(x) = Vq(x) — AV0(— x), where 
B(x) is the unit step function. The actual value of the step 
is much larger than the Fermi energy of the 2DEG. This is 
frequently used to set the potential step to infinity. However, 
here it is essential that the step is finite as the whole spin-orbit 
coupling is nonzero only in the region of the nonzero gradient 
of the confining potential. Since this value should be of the 
order of the work function of the electrons in the 2DEG, we 
fix this value to AV = 2.7eV— 230Ha*. One of the results of 
our work is the demonstration that the current induced spin 
polarization is rather independent of this value so that we do 
not need to be concerned with its precise numerical value, as 
long as it is large (>> £>) but finite value. 

Once we establish the analytical result for this model of 
abrupt potential step, we will also consider more general 
forms: step with a linear slope on a distance d, V c /(x) (used 
in Fig. [TJ, and a partially self-consistent, density-dependent 
model with a small triangular-shaped dip, V„(x), close to the 
edge. In the absence of this dip, the decrease in the density 
to zero at the edge in the 2DEG results in a slight local deple- 
tion of electrons there. A small negative value of the depth of 
the dip, —Vo, that is found self-consistently results in a charge 
neutral edge which is the physically expected situation. 

The essential part of our model is the spin-orbit (SO) cou- 
pling term in the electrons' Hamiltonian. In general, several 
different types of SO interactions can be similar in magni- 
tude: the Rashba-Bytchkov, the Dresselhaus or the impurity- 
induced SO coupling. However, since they are all small per- 
turbations to the effective Hamiltonian, we can consider them 
separably and superpose their outcomes in the sense of first 
order perturbation theory in terms of SO parameter appearing 
in the SO interaction. In our paper we will consider a special 



case of the impurity-induced SO coupling for the special case 
when the "impurity" is the edge potential only. In general, the 
SO contribution to the Hamiltonian takes the form 27 

F = a £ ff-(lxW(i)j, (2) 

where Cfe characterizes the strength of the SO coupling, a 
are Pauli matrices, k is the operator of electron's momen- 
tum and V(x) is the effective one-electron potential energy. 
Within our study the latter corresponds to various models 
of the edge potential: Vo(x), Vd{x) or V n (x). The strength 
of the SO coupling in GaAs is known to be approximately 
(Xe ~ 5.3A 2 = 5.53 x 10~ 4 ajj which is 10 6 times larger than 
the strength of the SO coupling in vacuum 2 ^. 

Within the framework of this model it is very easy to under- 
stand the origin of the current-induced spin polarization. The 
gradient of the potential energy is directed only in the x direc- 
tion, the momentum is confined only within the xy plane so 
that the only nonzero component of the spin operator comes 
from its z component due to the mixed product form in Eq. [2] 
Specifically, the SO interaction is simply a spin-dependent po- 
tential 

V so = -a E a z qV'{x) (3) 

with a well defined quantum number for the projection of the 
spin in the z direction, a z , attaining values ±1, and q is the y 
component of the electron's momentum (Fig. [TJ. Hence, the 
spin-orbit term is a potential energy which, for states with pos- 
itive momentum in the +y direction (the current), is smaller 
(larger) in the spatial region of negative derivative of the edge 
potential, V'(x) < 0, for electrons with their spin down (up). 
Since close to x = we have V'(x) = — AV8(x) < 0, we expect 
that close to the edge we will find majority of spin-down elec- 
trons. Of course, in equilibrium for zero total current density 
there will be equal number of electron with q > and q < 
so that in this special case zero total spin polarization will be 
observed. 



III. SPIN POLARIZATION AND THE SPIN-DEPENDENT 
PHASE-SHIFT 

In our previous work 2 ^, we have given estimates for the 
total induced spin per unit length of the edge using arguments 
based on wavepacket propagation. We considered a single- 
electron wavepacket state, V^'mM' f° r a specific momentum 
in the y direction q, spin a, with momentum k of uncertainty 

rk+Ak/2 dk / , , \ 

X( ,-/(» ! /2+,'-/2). i (4) 

where r q ' a (k) is the reflection amplitude of the eigenstate with 
energy E = (l/2)(k 2 +q 2 )- For large negative times only the 
incoming part (moving from the right in the Fig. [JJ of this 
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FIG. 1: (Color online) The edge's confinement potential of depth AV 
leads, via the SO coupling (Eq. [3j to attractive potential (V s0 in 
green) for electrons with spin down in the vicinity of the edge. The 
initially compensated spins of electrons approaching the edge be- 
comes locally un-compensated due to the difference in the scattering 
shift/ . 



of the edge, as given in our previous work. Here we will show 
that this physically motivated estimate is in fact a rigorous re- 
sult that we formally prove in the following, directly using the 
eigenstates of the Hamiltonian of the studied system. 

The spin polarization per unit length of the edge is given in 
terms of the spin-resolved density as 

m z = dxj^j ^(n\' q {x) -nj*(*)), (7) 

where the integration in q and k goes over all the occupied 
eigenstates and n ( f(x) is the contribution to the density at po- 
sition x.y from an occupied eigenstate 



„ — ikx 



(k)e 



= 2+29* 



j e /(2fa+tt«>") j 



(8) 

Since the spatial shift l qa is very small (given by the strength 
of the SO interaction), we can approximate the expression for 
the spin density to the first order 



wavepacket will give rise to nonzero contributions (the re- 
flected part will be extremely small due to the quickly oscillat- 
ing factors) and the wavepacket will be approaching the edge 
(Fig-E moving to the left). On the other hand, for large posi- 
tive times only the reflected wave will give nonzero contribu- 
tion. The position and form of the reflected wavepacket will 
depend on the spin of the incoming wavepacket since the re- 
flection amplitude, r qa {k) will be different for the two possi- 
ble spin orientations. Specifically, it is shown in the Appendix 
A that the reflection amplitude has the form r q ' a (k) - 
where the phase shift up to 3rd order in cce is 



q.a 



e. 



q-t/i _ 



2„2 



7z; + 2atan— ^2ka E q + 2kxa E q 
T 4a 3 E q 3 k(AV~^k 2 ), 



(5) 



where K — \J2AV — k 2 . Assuming the width Ak of the 
wavepacket Y^t a*vO sma H> one can easily find that the proba- 
bility density of the reflected wavepackets with the phase-shift 
given by Eq. |5]for large times will take a form 



1 sin 2 (x(t)Ak/2) 
% x(t) 2 Ak/2 



(6) 



where x(t) = x + d6 q,a /dk — kt, i.e. the maximum of the prob- 
ability is at x = kt — dQ c, ° jdk. From this it is evident that 
the electron with spin down will be lagging behind that with 
spin up by a distance 21" = d9 q ^ jdk — dO q '^/dk (indicated 
for the reflected right-going wavepackets in Fig. [TJ. This be- 
havior of the scattering of single electron in the wavepacket 
can be directly extended to many-electrons in view of the stro- 
boscopic wavepacket basis 22 . Namely, the above wavepack- 
ets for times t m = t + 2Km/(kAk), m = 0,±1,±2, where t 
is arbitrary moment of physical time, form orthogonal set of 
states which are all occupied with electrons. Hence the prod- 
uct of the length 21" times the number of electrons per area 
must give an estimate of the spin-polarization per unit length 



n ( f{x) = n c q (x) 



k.q 



i=o dl 



l q ' a , 



1=0 



n f(x) 



1 d 

+ 

1=0 2 dx 



nf(x) 



[>!■> 



(9) 
(10) 



where we have interchanged the differentiation in view of the 
functional dependence of the density, Eq. [8] on x and l q,a . 
Hence, using Eq.|7]and Eq.[l0]we obtain for the spin polariza- 
tion 



dq 
2k 



dk tg 

2% T 



(x) 



l q 

1 ,i ! 



(ID 



where we have introduced the antisymmetric part of the spa- 
tial shift 



0<7 t _ Qq,l 
2k 



= -2a E q-4alq 3 {AV--k z ). (12) 



The density at x = —°° is exponentially small, where as the 
contribution in the bulk of the 2DEG is according to Eq. [8] 
equal to 2. 

Including higher order terms in Eq. [10] would result in ap- 
pearance of contributions with higher spatial derivatives of 
the density at x = ±°°. There, however, is the density con- 
stant (either zero for x — —°° or the homogeneous 2DEG value 
for x = +°°) so that the expression for the spin polarization, 
Eq. QT|is correct also for large values of l q -° in principle. 

To complete the derivation we need to integrate over the 
occupied states. The simplest way is to assume that the non- 
equilibrium distribution is well described by a Fermi sphere 
with the radius (Fermi momentum) k E and the correspond- 
ing Fermi energy Ep — k F /2 in 2D, shifted by a drift mo- 
mentum qj in the +y direction. Such a model corresponds 
to the distribution function maximizing the information en- 
tropy for a fixed average energy, density and curren t 28 i 29 , and 
is also supported by Boltzmann-equation based analysis 30 . 
In the Appendix B we discuss the results obtained using the 
non-equilibrium model with two Fermi energies, conventional 
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within the meso- and nano-scopic quantum transport. It gives 
identical result for the first order expansion in terms of a E ; 
higher order terms in a E differ by a numerical prefactor but 
are negligibly small for any realistic situations. 

Using the expression for the asymmetric spatial shift, 
Eq. [T2l we find for the induced spin polarization linear in the 
drift momentum (and hence the current density) 

= (-a E -3a E AVE F + ^a E E F ^q df i 2D (13) 

where kp — y/2Ep is the Fermi momentum, «2D = Ep/n is 
the density of the 2DEG. (See Fig.|5]for the explanation of the 
integration domain for k, q.) On the other hand, the current 
density is j = q ( in2D so that we find the final result 

m z = -a E j -3a E AVE P j + ^a E E F j + &(a E ,f) (14) 

In a typical situation, all except the first term of this expan- 
sion are negligibly small so that the resulting spin polarization 
is independent of the height of the confining potential AV and, 
according to the numerical results presented in the following 
section, practically independent of other gentle modifications 
of the model so that the first term should be useful as a general 
estimate of the magnitude of the current induced spin polariza- 
tion in heterostructures based on 2DEG. Expressing the first 
order result in the S.I. units we obtain the central result of our 
work, stated in the introduction in Eq. [T] 

IV. EFFECTS OF THE FINITE SLOPE OF THE EDGE 
AND THE PARTIAL SELF-CONSISTENCY 

The aim of this section is to explore the rigidity of the result 
in Eq. Q]with respect to changes in the model of the confining 
potential at the edge. We will first address the dependence of 
our result on finite slope of the edge potential on the distance 
d using the potential energy 

( AV x<0 

V d (x)=l AV-(AV/d)x 0<x<d (15) 
[ x > d 

Qualitatively we can expect this dependence to be similarly 
weak as the dependence on AV obtained in the previous sec- 
tion (Eq.[T4"l). The reason behind this is a mutual cancellation 
of two effects: (1) increasing AV increases the strength of the 
SO coupling in the Hamiltonian and hence the spin polariza- 
tion, (2) increasing AV decreases the amplitude of the wave- 
functions in the region of the nonzero gradient of the confining 
edge potential and hence its sensitivity to the SO coupling. 

To characterize this dependence quantitatively we have 
evaluated the expression Eq. |7] numerically, where the spin- 
dependent contribution to the density, n a q (x\ was obtained 
from the eigenfunctions of the effective ID Hamiltonian con- 
taining the potential energy Eq. Q3] and, according to Eq. [3] 




FIG. 2: (Color online) The difference in spin densities (red-full) 
and the density (green-dashed) of electrons close to the edge of the 
sample (located at x = 0). The former, similarly to the density, ex- 
hibits Friedel-like oscillations. Its first dominant local minimum (at 
x ~ 0.2) dominates the contribution to the nonzero spin-polarization 
obtained by integrating the difference in spin densities over whole 
x axis. The inset demonstrates the extrapolation of this integration 
with respect to its upper limit, /, using a fit to a functional form 
~Asin(2k F l + <p)/l. 



from it derived spin-orbit potential energy. The eigenfunc- 
tions were obtained by matching the plane-waves in regions 
x < and x > d with the Airy functions in the region of the 
potential energy with linear slope, < x < d. The integration 
over k in Eq. [7] could be done analytically, whereas the inte- 
gral over q was done numerically on the regular mesh with 
N/i points. We used model parameters that are typical for 
2DEG created within GaAs quantum wells?, electronic den- 
sity n = 0.985 corresponding to the Fermi energy Ep =3.01, 
the current density j — 0.1 and the value of the spin-orbit pa- 
rameter a E — 5.53 x 10~ 4 . 

For if = 0we confirm our analytical results for abrupt step 
potential, Eq. [14] The difference between the spin-up and 
spin-down densities, shown in Fig. |2] exhibits Friedel-like 
oscillations with the first period having a pronounced nega- 
tive amplitude giving the major contribution to the non-zero 
spin polarization. The proportion of this amplitude to the 
electron is ~ 10~ 4 for the current density j = 0.1. Since 
the current density in the experimental situation is typically 
j = 0.1 — 0.01 (the amplitude scales linearly with the current 
density for these values of j), this value is comparable to the 
amplitude of the spin polarization in the extrinsic spin Hall 
effect 4 . Unfortunately, in contract to the extrinsic spin Hall ef- 
fect, the nonzero spin-density is located only within few Fermi 
wavelengths away from the edge. 

Integrating the difference in the spin densities along the 
whole x axis we obtain the spin-polarization per unit length of 
the sample. The numerical result of this integration naturally 
depends on the upper limit of the integration; a reliable result, 
is easily obtained extrapolating the upper limit to infinity (in- 
set of Fig. [2]). Using this methodology we find the numerical 
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FIG. 3: (Color online) The dependence of the self-consistent value 
of the potential dip Vq on the electronic density shows monotonic be- 
havior. While the value of the dip is comparable to the Fermi energy 
(Ef = 3.01), the considered values do not lead to appearance of the 
bound state at the edge. 
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result for the spin polarization m z = (5.53 ±0.003) x 10 5 , 
which is in perfect agreement with our analytical expression, 
Eq.M 

Using this extrapolation scheme we can address the 
changes in the spin-polarization with parameters of the model 
discussed in the following. First, we confirm the observation 
that the dependence on the magnitude of the confining edge 
potential comes in the 3rd order of Ofe: the value of the slope 
of the polarization vs. the height of the confinement edge po- 
tential, AV is dm z /d(AV) 1.6 x 10~ 10 which is close to the 
analytical value 1.5 x 10~ 10 . The difference comes primarily 
from the extrapolation of the value of m z with respect to the 
upper limit of integration in x. 

Considering nonzero values of the smearing length of the 
confining potential, d € (0,Xf/2), we find that the resulting 
spin polarization is changing only very little (dm z /dd ~ 3.3 x 
10~ 8 ). The reason for this negligible dependence is, similarly 
to the dependence on AV, the mutual cancellation of the two 
effects discussed earlier. 

The second modification to the confining potential that we 
consider is partial self-consistency of the edge potential that 
guarantees charge neutrality of the edge of the sample. Pres- 
ence of the model confining edge (Eq. IT3T > leads to redistribu- 
tion of the charge density characterized with excess positive 
charge (per unit length of the sample) 

AN = j + dxAn(x) = J + dx{E F /%-n{x)), (16) 

where n(x) is the electronic density and Ef jn is the density 
of the positive background charge. The Fermi energy of the 
electrons guarantees that for large x we have n(x) = Ef /n so 
that the system is charge neutral in the bulk even through the 
total charge at the edge is not necessarily zero. The behav- 
ior of the density at the edge of a 2D sample is in contrast 
with the typical situation for the edge of 3D metals where the 
work function and the Fermi energy, mutually comparable in 
value, lead to leakage of the electronic density into vacuum 
and thereby to an access of negative charge 31 . However at the 
edge of a 2DEG, the Fermi energy, typically few meV is much 
smaller than the work function which is of order of several eV 
and the situation is reversed. 



FIG. 4: (Color online) The difference in the spin densities for non- 
interacting (red-full) and the partially self-consistent (green-dashed) 
electrons close to the edge of the sample (located at x = 0). The 
spatially resolved spin density difference exhibits pronounced dif- 
ferences due to the selfconsistency, the most important effect being 
the shift of the curve closer to the edge. However, this results only 
in weak enhancement of the total spin polarization, m z , shown in 
the inset for an interval of electronic densities, due to counter-acting 
spin-orbit coupling induced by the self-consistency correction in the 
region x <E (0, Xf /2). 



The overall charge of the edge can be neutralized by a sim- 
ple model potential 

( AV x<0 

V„(x) = \ -V + (2Vo/X F )x 0<x< Xf/2 , (17) 
[ x>X F /2 

which on a distance Xf/2 exhibits small potential dip to attract 
more electrons. The magnitude of the dip is obtained self- 
consistently to achieve zero total charge at the edge, i.e. the 
condition AN = 0. While in principle this form of the potential 
may support new bound states (edge states) for sufficiently 
large Vq, we have checked that for the here, self-consistently 
found values of Vq no such states exist. The dependence of the 
self-consistent value of Vq on the typical values of the density 
of the 2DEG is shown in Fig. [3] 

Including the SO interaction, Eq. [3] on top of the potential 
energy V„(x), we obtain aditional repulsive SO-induced con- 
tribution for the spin-down electrons due to the linearly rising 
potential of the dip for < x < Xf/2. This is in competition 
with the character of the edge-confinement potential prefer- 
ring the spin-down electrons. However, the dip in the potential 
energy also tends to move the electrons closer to the confining 
edge and, by increasing the density in this region, effectively 
enhances the effect of the SO coupling there. While the form 
of the spin-density does change significantly due to all these 
mechanisms (see Fig. [4j, the overall spin polarization per unit 
length remains essentially unaltered over a wide range of elec- 
tronic densities, shown in the inset of Fig. [4] This result once 
again confirms the rigidity of the result given by Eq. Q] 
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V. CONCLUSIONS 

The confining potential close to the edge of a 2D electron 
gas together with a nonzero current along this edge induces a 
nonzero spin polarization that is localized within a few Fermi 
wavelengths from the edge. To characterize this effect quan- 
titatively we have derived a simple analytical formula for the 
spin polarization per unit length of the edge. Interestingly, the 
spin polarization is independent of the height of the confining 
potential as well as the electronic density of the 2DEG. Fur- 
thermore, using numerical calculations we have also showed 
that this result is independent on other possible modifications 
of the shape of the confining potential: the spatial extent of the 
confining potential and the partial selfconsistency of the con- 
fining potential with respect to charge neutrality of the edge. 
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FIG. 5: The shifted Fermi sphere (left) and the two-Fermi radii 
(right) occupations in 2DEG. The occupied states are shown as filled 
areas, have both identical current density and electronic density but 
differ in higher moments of k. 



Expanding the last result in powers of a E to 3 rd oder we obtain 
the expression Eq. [5] 
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VI. APPENDIX A 

The calculation of the spin-dependent phase-shift follows 
the usual textbook treatment of the wave-function matching 
method in ID problems. For the separated x~ dependent fac- 
tor of the eigenstate we need to solve the ID Schrodinger 
equation 



AV9{-x)±a E qAV8(x) $f a {x) = e^ a (x), 



Id 2 

'lib? 

(18) 

where the '+' sign is for a =f spin and the ' — ' sign for a =\. 
spin. We are only interested in the solution well below the 
vacuum, i.e. e < AV for which we demand asymptotic forms 



^■ a (x)=t q - a {K)e KX , x<0,K= v / 2{AV-e) (19) 
$>°(x) = e - ikx + r q - a {k)e ikx , x>0,k=V2^ (20) 

The reflection amplitude (r q a (k)) and the coefficient t q ' a {K) 
are then found from the continuity of <j>e ' (x) and its derivative 
at x = with the result 



r q - a {k) = - 
t q ' a (K) = - 



±2a E qAV + K + ik 
±2(XEqAV + K — ik 
2ik 

±2a E qAV + K + ik' 



(21) 
(22) 



The sought phase shift is obtained from the reflection ampli- 
tude rf a = \r4' a (k)\e i0 k a , 



VII. APPENDIX B 

We have stated that the result for the spin-polarization, 
Eq. [13 is independent of the considered non-equilibrium oc- 
cupations of electronic states to the first order in the the spin- 
orbit coupling <Xe- This is true as long as the current density 
in the system is the same for all of these considered occupa- 
tions. Let us consider two simple models of occupations: (1) 
the "two-Fermi radii model" (2FR), frequently used within the 
coherent transport with the occupancies dictated by the left 
and right macroscopic electrodes with Fermi momenta k F and 
k F respectively 32 attached to the sample at its ends (Fig. [5] 
right), and (2) the shifted Fermi distribution function (ShF), 
used within the main text (Fig. [5] left) From this one immedi- 
ately sees that the contribution to the spin polarization that is 
linear in a E in Eq.Ql] and therefore also linear in the momen- 
tum k y (— q), will be independent of the particular form of the 
occupations as long as the zero-th and the first moments are 
identical. 

Apart from this, the 2FR and ShF occupations represent two 
different situations: the former is suited for very short ballistic 
systems and it directly facilitates interpretation of the results 
in terms of the difference in electro-chemical potentials of the 
electrodes, Afi = (1 /2) ((kp) 2 - (kf) 2 ). The latter is related 
to the current density inside the sample through 



(24) 
(25) 



9 qa = 7z; + 2atan . 

* ±2a £ ^AV + k 



(23) 



where the occupation factor, n(k x ,k y ) = 2 (spin degeneracy 
in the unperturbed system) for the occupied states shown in 
Fig.|5]and n(k x ,k y ) = otherwise. Eq.|25]gives the linear ex- 
pansion in the applied bias Aji. No such simple connection to 
the applied bias voltage can be made for the ShF occupation. 
On the other hand, ShF is suitable for longer samples, where 
electrons' scattering leads to partial equilibration of the dis- 
tribution function 30 . This is then characterized with the drift 



7 



momentum, q c /, which gives the current density 



j = n 2 Dqd = — qd- 



(26) 



Since the samples we consider are typically longer than the 
electron's coherence length, we have used the ShF occupa- 
tions within the main text. For completeness we give also the 
results for the 2FR case. The spin polarization per unit length 
of the sample to third order in Ofe is 



2alAV 

— / 

3n 2 



,3/2 



am 



4a 



45n 2 



(27) 



which after expressing the applied bias in terms of the current 
density results in 



-«EJ- 



2alAV , 



-E F j - 



4a 



E p2 



45 



fJ- 



(28) 



Comparing the last equation with Eq.[l4]we see that while the 
first order term is identical for both occupations, the higher 
orders differ by a numerical prefactor. Both of these are larger 
in the case of partially equilibrated shifted Fermi-like occupa- 
tions. 
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